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Abstract 

The hole distribution in a double quantum well (QW) structure consisting of a magnetic and 
a nonmagnetic semiconductor QW is investigated as a function of temperature, the energy shift 
between the QWs, and other relevant parameters. When the itinerant holes mediate the ferromag- 
netic ordering, it is shown that a bistable state can be formed through hole redistribution, resulting 
in a significant change in the properties of the constituting magnetic QW (i.e., the paramagnetic- 
ferromagnetic transition). The model calculation also indicates a large window in the system 
parameter space where the bistability is possible. Hence, this structure could form the basis of a 
stable memory element that may be scaled down to a few hole regime. 

PACS numbers: 72.20.Ht,85.60.Dw,42.65.Pc,78.66.-w 
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Spintronicsi approaches the brink of apphcations in the micro- and nano-electronics. As 
soon as the semiconductor compounds doped with magnetic ions reveal room temperature 
ferromagnetism, the creation of competitive electronic devices will become the barest neces- 
sity. The study of different structures and their spin effects, which could result in spin devices 
at lower temperatures, is already a topical problem especially when we take into the account 
the continued improvement in the critical temperature Tc of the paramagnetic-ferromagnetic 
(PM-FM) phase transition. 

In this Letter, we explore the formation of a bistable state and its device application in 
a double quantum well (QW) semiconductor structure with respect to the distribution of 
itinerant holes between the constituting nonmagnetic and magnetic QWs. The calculation 
illustrates that the interplay between the free carriers and the magnetic ion spins is the key 
to achieving the bistability where the magnetic QW can switch between the PM and FM 
phases in a controlled manner. 

The structure under investigation is illustrated in Fig. 1. Consider a magnetic QW 
(MQW) and a nonmagnetic QW (NQW) of widths L^m and L^n, respectively, separated 
by a barrier. The width Lb and the hight of the barrier should be large enough to form 
non-coherent single QW states and yet not too large in order to enable hole redistribution 
when a gate bias is applied. The total 2D hole concentration n° is assumed to be a constant 
and the FM-PM transition in the MQW to be mediated by two-dimensional (2D) free holes. 
If the hole energy in the NQW is lower than that in the MQW, one can expect that the 
stable state will correspond to hole localization primarily in the NQW with a small leakage 
in the MQW, which is in a PM phase (see the schematic on the left in Fig. 1). When a 
proper bias is applied, holes from the NQW can be transferred to the MQW via tunneling, 
over-barrier injection, etc. As the hole density in the MQW surpasses a certain threshold at 
a given operating temperature, the layer undergoes the PM-FM transition. When the hole 
exchange interaction with the FM-ordered ion spins in the MQW is strong, it can reduce 
the total free energy below that of the initial state with the PM phase even after the bias 
is switched off. Then, the holes will remain confined in the MQW and the FM state be 
maintained (the schematic on the right in Fig. 1). When a reverse bias pulse is applied, 
the holes are drained out of the MQW into the NQW, and the MQW will return to the 
PM state. Hence, if realized, these two stable states can coexist under the same external 
condition and the switching between them is mediated by the electrically-controlled PM-FM 
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phase transition.— It is expected that the structure can operate up to a temperature shghtly 
below the saturated maximum of Tc, which can reach room temperature or higher in some 
material systems.— 

To analyze this problem, typically one would derive the magnetic Hamiltonian H„i, which 
describes the effect of the effective ferromagnetic inter-ion spin-spin interaction in the pres- 
ence of free holes. Then the calculation of the free energy F of the total system consisting 
of magnetic ions and free carriers, which occupy the MQW with a concentration of rihM and 
the NQW with n/^Tv iP'hM + = n^h)i leads to the carrier population factor 77 = nhM/n^ 
and 1 — 1] = UhN/n^ at each QW that minimizes F = F{r]). The existence of two local 
minima with respect to rj demonstrates the bistability of the system under consideration. 

Although conceptually correct, this approach faces the difficulty of specifying the mecha- 
nisms responsible for FM ordering in the MQW (for details on various mechanisms, refer to 
those cited in Ref. To circumvent this problem, which is beyond the scope of the present 
study, we develop a semi-phenomenological approach that utilizes the data extracted from 
routine experimental measurements of magnetism. Namely, we assume that the magnetic 
part of the free energy Fm can be expanded with the magnetization M (order parame- 
ter) according to Landau theory. This expansion approximates satisfactorily the veritable 
dependence in the whole temperature interval we are interested in. Assuming the easy 
magnetization axis is directed along the growth axis of the sample and a magnetic field ^ 
parallel to M, the Fm expansion in the most general form reads 

Fm = -a{T, - T)M^ + hM^ - ^MB. (1) 

It can be shown that the parameters a and b of Landau expansion in Eq. (0) can be expressed 
in terms of fundamental properties of the magnet: The Curie- Weiss law for the magnetic 
susceptibility x = Cq/{T — at T > Tc defines a = I/ACq, while the spontaneous 
magnetization Mg = [a(Tc — T)/2bY^'^ = MqJi — T/T^ for T < (that minimizes Fm 



at 5 = 0) provides b = aTc/2MQ. By these relations, all parameters in Fm in Eq. 
are determined. Tc can be found separately as a function of the magnetic layer carrier 
population. 

Since we are looking for the free hole distribution over the magnetic and nonmagnetic 
QWs at i? = 0, the magnetization M in Eq. (P) takes the equilibrium value. Hence, 

a\Tc-Tf 
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T < Tc, Fm = 0, T > Tc. (2) 



It is important to note that Eq. (j2I) includes the dependence on the free hole concentration 
in the MQW via the critical temperature Tc = Tc{ri). 

The total free energy of the system can be obtained if Eq. ^ is supplemented with the 
free energy of the hole gas: 

F,, = FMv) + U{v)+C{v). (3) 

The first term in Eq. ^ accounts for the kinetic energy of the hole gas in both QWs, U{r]) 
is the hole potential that is different for the MQW and the NQW, and C{ri) is the energy 
of the Coulomb interaction between the QWs. Assuming the parabolic dispersion law with 
an effective mass m for 2D holes and that only the lowest sub-bands of the QWs can be 
populated by holes with L^m = L^n, one can find the kinetic energy for the hole gas as 

F2D{ri) = ksT jr/A (^^r^^ + (1 - r^)^ (^^(1 - r^)^ | (4) 

where ks is the Boltzmann constant, E^p = T:h^n\/m is the Fermi energy of 2D holes with 
concentration ?t,°, and 

= In (e^-l) + 1^^2(1 -e^) (5) 

X 

with a polylogarithmic function Li2 (x) = /°(itln(l — t)/t. At low temperature T Sp/k, 
Eq. (j3)) describes the sum of degenerate carrier energies in both QWs. The low temperature 
assumption is commonly used in the works on the FM ordering in MQWs. For the present 
purpose, however, we need to account for the arbitrary relation between e% and /c^T accord- 
ing to Eq. (0]). The contribution of the energy shift AU between the MQW and the NQW 
(Fig. 1) is accounted for in the term 

U{r]) = Af/r/<, (6) 

while the Coulomb energy in the strong confinement limit takes the form 

CM = ?^L..r(,-i)\ (7) 

where e is the electron charge and e the dielectric constant. 

Now one can analyze the total free energy F = Fm + F^ with respect to possible FM phase 
transitions in the MQW. Considering that the realistic dependence Tc = Tdrihu) should be 
taken from the experiments, we specify our analysis by a model that can be applied to 
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a typical diluted magnetic semiconductor. Assuming that the total magnetization of the 
MQW stems mainly from the magnetic ions and n\ <^ rimLy^M {nm is the 3D magnetic ion 
concentration in the MQW), one can easily find the parameters a = 3/[4:S{S + l)g'^^\n^ 
and h = 3kBTc/[8S^{S + l)g'^fi%n^], where g is the magnetic ion (jf-factor with spin S and 
fiB denotes the Bohr magneton. 

In order to describe the dependence of Tc on rj, we propose an approximation 

T, = T° (l - e-4WA.BT0^ ^ 

where T° is the asymptotic (at a high enough UhAi) value of the critical temperature and 
a is the fitting parameter that adjust the dependence [Eq. (jH))] to the experiments. In 
the following, we assume a = 1 since it describes the experimental results satisfactorily.— 
Combining this approximation with Eqs. Q, and ((7j), we find the trial function in 
the form of the free energy per hole normalized by the energy unit fc^T^ 

F--^^^[te(.)-tf^(t.(,)-t) 

+U7] + wi27] - ly + F,^, (9) 
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where u = rimL^/nl t = T/T^, t,{r]) = TJT^, u = AU/knT^, w = TTe^Lbnl/2ekBT^ 
e^p/ksT^i tc{i]) = 1 — exp(— rr/), and 6{x) is the Heaviside step function. For completeness, 
Eq. © also includes the exchange energy of free carriers Ff,^, which can influence the effect of 
bistability.i^ We performed the calculation of F^x in a Hartree-Fock approximation following 
Ref. y, where the exchange and the correlation potentials of a 2D gas were found at finite 
temperatures. The final result of our calculation has a proper analytical approximation in 
terms of the dimensionless parameter r 

25/2e%/^ 
F — 1 - X 

v'^'4>{^v) + ii-v?^'<p{jii-v))], (10) 

where (x) = 1 - e'^^^-^ + xV(2.4 + x"), 6 = 0.3 + 0.29{x - 1), c = 2.2 - 0.2^(a; - 1). 
Note that Eq. takes a form similar to the free energy expression used for analyzing the 
spin/charge separation of magnetic semiconductors near a PM-FM phase transition.— 
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For a numerical evaluation, let us assume the following "typical" values for the parameters 
of the double QW structure: m = 0.3mo (mo is the free electron mass), e = 12.9, 5* = 5/2, 
L^M = L^N = 10 nm, U = 5 nm, nJJ = 10^^ ^,^^-2^ ^ ^ 3 ^ iq21 ^^-3^ ^^^^ rpo ^ ^qq 

Figure 2 displays F{r]) calculated at three different values of the energy shift AU between 
the minima of PM QW and NQW (m = AU/keT^ = 29,5,-9). Clearly, curves 1 and 3 
{u = 29, —9) support only one stable state at 77 = or 1. This means that when the MQW 
state lies either too high (curve 1) or too low (curve 3) compared to that of the NQW, the 
holes strongly prefer to be confined in one of the QWs. Even when they are transferred to 
the other QW through an external bias, the holes will return to the preferred state once the 
bias is turned off. However, one can realize a structure that has free energy minima at or 
near both rj = (with the MQW in the PM phase) and rj = 1 (the MQW in the FM phase) 
if AU is properly selected (curve 2). These two states can be stable with respect to small 
fluctuations under the same external conditions. We also found that the relative effect of 
Fex is small compared to the Coulomb energy of the inter-QW carrier interaction w{2ri — iy. 

This bistability can be achieved in a relatively wide range of AU and T as shown in 
Fig. 3. As expected, the condition for AU becomes less stringent with a decreasing T. This 
is because the structure can now operate with a lower Tc, which in turn requires a smaller 
hole density in the MQW for the PM-FM transition. The highest operating temperature 
will be somewhat lower than T° and is a function of the maximum possible n^M- 

It should be emphasized again that the structure is analyzed with respect to the free 
hole distribution r) that minimizes the system free energy. Hence, the height of the local 
maximum in the free energy (e.g., near rj = 0.8 for curve 2 in Fig. 2) does not constitute the 
energy barrier separating the NQW and the MQW. In fact, details of the barrier layer shown 
in Fig. 1 is not considered in the present study. Rather, the stability of local minima (thus, 
the lifetime) depends explicitly on the magnitude of the fluctuation in 77 that can cause 
unwanted switching. Since the required value is large according to our calculation (e.g., 
Ar] > 0.2), the system is expected to be robust against most fluctuations including thermal 
transitions. Hence, the lifetime at each local minima may be as long as that of conventional 
magnetic memory cells. A detailed analysis is necessary for a quantitative estimate. 

In summary, we demonstrate bistability formation in the structure consisting of magnetic 
and nonmagnetic QWs. The bistability is mediated by the PM-FM transition in the MQW. 
In contrast to the case of coupled NQWs,-^ our structure reveals a bistability as soon as the 
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conditions for the ferromagnetism in the MQW are satisfied. As a result, a room temperature 
operation may be achieved once a proper magnetic material is developed. Although the 
investigation is done for a 2D structure, it is expected that a similar principle can also be 
applied to the OD system. As the size of the gate electrode shrinks, the MQW can form a 
gated quantum dot, which can undergo the PM-FM (or super-PM) transition by controlling 
population/depopulation of holes. Furthermore, one can envision a magnetic nanocrystal 
embedded in a nonmagnetic barrier in place of the MQW. Hence, this structure can form 
the basis of a stable memory element that may be scaled down to a few hole regime with 
very low power consumption. 

This work was supported in part by the DARPA/ONR and the SRC/MARCO Center on 
FENA. 
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Figure captions 



Fig. 1. Schematic energy diagram (valence band) of the structure in the two coexisting 
stable states. The mutual ahgnment of the hole spins (large arrows) and locahzed spins 
(small arrows) reduces the energy in the FM QW by £'exch — SF^/n^ due to their exchange 
interaction. The switching between these states (arch arrows) is achieved by applying an 
appropriate bias pulse on the gate electrode. 

Fig. 2. Free energy trial function F{r]) at different values of u (= AU/ksT^)- Three different 
scenarios are shown: curve 1 {u — 29) - a monostable case with holes occupying the NQW 
(77 = 0), while the MQW is in a PM phase; curve 2 {u — 5) - a, bistable case where the PM 
and FM phases coexist; curve 3 {u — —9) - a monostable case with holes populating the 
MQW (77 = 1) in the FM phase. 

Fig. 3. Phase diagram of the parameter space indicating the potential bistability region. A 
free hole density of 10^^ cm~^ is assumed. 
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Fig. 1: Semenov et al. 
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Fig. 2: Semenov et al. 
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